Topological phases in layered pyrochlore oxide thin films along the [111] direction 
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We theoretically study a multi-band Hubbard model of pyrochlore oxides of the form A2B2O7, 
where B is a heavy transition metal ion with strong spin-orbit coupling, in a thin film geometry 
orientated along the [111] direction. Along this direction, the pyrochlore lattice consists of alternat- 
ing kagome and triangular lattice planes of B ions. We consider a single kagome layer, a bilayer, 
and the two different trilayers. As a function of the strength of the spin-orbit coupling, the direct 
and indirect ci-orbital hopping, and the band filling, we identify a number of scenarios where a 
non-interacting time-reversal invariant Z2 topological phase is expected and we suggest some can- 
didate materials. We study the interactions in the half-filled d-shell within Hatree-Fock theory and 
identify parameter regimes where a zero magnetic field Chern insulator with Chern number ±1 can 
be found. The most promising geometries for topological phases appear to be the bilayer which sup- 
ports both a Z2 topological insulator and a Chern insulator, and the triangular-kagome-triangular 
trilayer which supports a relatively robust Chern insulator phase. 



PACS numbers: 71.10.Fd,71.10.Pm,73.20.-r,73.43.-f 
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I. INTRODUCTION 

With the experimental discovery of two 1 ^ and three-- 
dimensional topological insulators, the study of topo- 
logical phases has rapidly intensified over the past few 
years 1^— Significant effort has been directed at the iden- 
tification of new materials^— and the role that interac- 
tions may play in driving topological phasesJ^— In ad- 
dition, the interplay between magnetism and non-trivial 
topological bands has also been addressed^ - — 

A number of the studies focused on interaction ef- 
fects have been directed at topological phases in transi- 
tion metal oxides £2r— A new direction that has emerged 
within this area is the search for topological phases at 
the interfaces of correlated oxides 42r— 

To date, the prediction of topological phases in tran- 
sition metal oxide interfaces has focused on the per- 
ovskite structure 4 ^— ABO3, where A is usually a rare 
earth element, B is a transition metal, and O is oxygen. 
An undistorted perovskite has a relatively simple cubic 
structure with natural cleave planes along the [001] and 
equivalent directions. For this reason, most of the inter- 
face and supcrlattice structures of these materials have 
been grown along this direction^ - — However, theoreti- 
cal models suggest that thin films, particularly bilayers 
and trilayers, grown along the [111] direction are more 
promising for realizing topological phases42r— The un- 
derlying reason for the increased favorability of topologi- 
cal phases in thin films grown along the [111] direction is 
the nature of the energy bands (derived primarily from 
the transition metal d-orbitals and the oxygen p-orbitals) : 
They are both relatively "flat" (i.e. dispersionless) along 
some directions for certain filling fractions, and they pos- 
sess flat band touching points in relevant tight-binding 
models. In two dimensions, flat band touchings are 
perturbatively unstable to interactions and topological 
phases often emerge as the leading instability] 23 : 50 : 52 : 63 : 64 



Flat bands themselves harbor topological phases if their 
Chern number is non-zero, including the most interesting 
case of partial band filling^ 5 - - — 
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(a) Sandwich 
structure 



(b) Brillouin zone 



(c) Pyrochlore lattice structure 



FIG. 1. (color online) (a) Schematic of the experimental ge- 
ometry. The middle A2B2O7 thin-film layer is the region of 
physical interest. The capping regions A2B'2 07 would ide- 
ally be non-magnetic large gap band insulators of the same 
lattice constant as A2B2O7 and would only serve to struc- 
turally stabilize the thin-film region. The growth direction 
of the structure is along [111], (b) The Brillouin zone of the 
underlying triangular Bravais lattice of the thin-films, (c) 
Pyrochlore lattice structure showing alternating kagome and 
triangular lattice planes along the [111] direction. A kagome- 
triangular-kagome-triangular tetralayer is shown with green 
atoms in the kagome planes and gray atoms in the triangular 
planes. 



A potential experimental challenge in the implemen- 
tation of the theoretical proposals^— for topological 
phases in the perovskite ABO3 heterostructures grown 
along [111] is the difficulty of crystal growth in this di- 
rection. As it is not a natural cleavage plane, it is 
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also not a natural growth direction. It is unclear at 
present whether this issue will become a significant im- 
pediment to the exploration of topological phases in ox- 
ide thin films and heterostructures, although successful 
growth of (lll)-oriented ABO3 heterostructures was re- 
cently reported.— ~— In the meantime, it is prudent to 
suggest other material systems where strong correlation 
effects and topological phases are likely to coexist. 

In this work we study the potential for realizing topo- 
logical phases in thin films [see Fig. [la] of transition metal 
oxides of the form A2B2O7, where A is a rare earth ele- 
ment, B is a heavy transition metal ion with strong spin- 
orbit coupling, and O is oxygen. In this class of materials, 
the energy bands are expected to be derived primarily 
from the transition metal ion d-orbitals and the oxygen 
p-orbitals ) 40 ' 42 ' 43 ' 48 similar to the ABO3 perovskites dis- 
cussed just above. However, in the case of A2B2O7 the 
transition metal ions B reside on a pyrochlore lattice-a 
corner sharing network of tetrahedra shown in Fig. Qx] 
Along the [111] and equivalent directions there are nat- 
ural cleavage planes that consists of alternating layers of 
a triangular lattice of B ions and a kagome lattice of B 
ions. [Sec Fig. [TcJ] Therefore, in contrast to the ABO3 
materials, the [111] direction is expected to be a good 
growth direction. 




(a) Kagome layer (b) Bilaycr 




(c) Trilayer (TKT) (d) Trilayer (KTK) 



FIG. 2. (color online) Candidate thin films to be used in the 
A2B2O7 region shown in Fig llal As shown in Fig llcl along 
the [111] direction the pyrochlore lattice is an alternating se- 
ries of kagome and triangular lattice planes. Single, double, 
and triple layer stackings are shown in (a)-(d). Among these 
stackings, the bilayer and TKT trilayer are found to be the 
most promising for realizing robust (to model Hamiltonian pa- 
rameters) topological insulator and Chern insulator (zero 
magnetic field quantum Hall) phases. 

In this paper we explore the favorability of topological 
phases in heterostructures grown along the [111] direc- 
tion with a geometry such as that shown in Fig. [Ta) The 
"bread slices" of the "sandwich" composed of the ma- 



terial A2B'2C>7, would ideally be a non-magnetic band 
insulator which would differ from the thin-film layer only 
in the "B" element to minimize lattice strain effects.— 
We consider a central thin-film layer that could consist of: 
(i) a single kagome layer, (ii) a triangular-kagome bilayer, 

(iii) a triangular-kagome-triangular trilayer (TKT), and 

(iv) a kagome-triangular-kagome trilayer (KTK). These 
thin film structures are shown in Fig. [2J We study a 
multi-band Hubbard model with an on-site atomic spin- 
orbit coupling of the L ■ S form, Eq. (pj. The inter- 
actions are included with the standard on-site Hubbard 
term, Eq. ((3]), and we treat them within the unrestricted 
Hatrcc-Fock approximation. 

Our main results are that we find <i-shell filling frac- 
tions in all systems studied that favor Z 2 time-reversal 
invariant topological insulators (TI) at small and van- 
ishing interactions. The specific filling fractions in each 
system that support a TI phase arc shown in Figs. [3][6j 
along with representative band dispersions along high- 
symmetry directions for two values of spin-orbit coupling 
and two values of hopping parameters. In order to com- 
pare our results for the case of arbitrary interactions with 
those obtained in bulk A2B2O7 pyrochlore a 46 ' 48 we focus 
on the case of strong spin-orbit coupling and a half-filled 
rf-shcll. Our main results for the non-interacting case are 
summarized in the phase diagrams in Fig. [7] and for the 
interacting case are summarized in the phase diagrams 
in Fig. [8j Two results are worth highlighting: (i) The bi- 
layer possesses a TI at small interactions and a narrow re- 
gion of a zero magnetic field Chern insulator with Chern 
number ±1 at intermediate interaction strength, (ii) The 
TKT system possesses a relative wide range interactions 
of intermediate strength where a Chern insulator phase 
with Chern number ±1 is stabilized. 

In our study, we have swept out a wide portion of 
the parameter space accessible within the tight-binding 
model and Hartree-Fock calculations employed (espe- 
cially in the case of a half- filled d-shcll). It is impor- 
tant to list some materials that might serve as a useful 
starting point in an experimental search for the topolog- 
ical phases we find here. For half-filled d-shell materi- 
als, materials such as A2h"207 with A=La, Y might be 
good candidates for the thin films considered here. For 
the "capping" layers A2B'207, possible candidates are 
A 2 Hf 2 7 for A=La, Y^I 8 . Undoubtedly there are oth- 
ers, and we hope material growth experts will pick up 
the thread and apply their own expertise to suggest the 
most promising candidates ! 79 ' 80 

Our paper is organized as follows. In Sec. |H] we 
describe the effective tight-binding model we consider, 
which includes an atomic spin-orbit coupling term that 
drives topological insulator phases and both direct d- 
orbital hopping between transition metal ions and the 
more commonly considered indirect d-orbital hopping via 
the oxygen p-orbitals. We present phase diagrams and 
band structures for a variety of representative thin films. 
In Sec. IIIII we include an additional on-site Hubbard in- 
teraction that we study within the Hartrcc-Fock approx- 



3 



imation. Wc present phase diagrams for the case of half- 
filled d shells for a variety different thin films and high- 
light the conditions that favor interaction-driven topo- 
logical phases with a non-zero Chern number. In Sec. IIVI 
we present the main conclusions of our work, and in the 
appendicies we provide some technical details related to 
computing the direct d-orbital hopping parameters on 
the pyrochlorc lattice (including the thin film geometries 
we consider) and the mean-field calculation. 



II. EFFECTIVE TIGHT-BINDING MODEL 

As we mentioned in the introduction, we are focused 
on the physics of the thin film A2B2O7 layer shown in 
Fig. [la] We consider the four varieties of films shown 
in Fig. [5] The "capping" layers A2B'2C>7 are assumed 
to be electronically inert, and this assumption has been 
substantiated in closely related density functional theory 
(DFT) studies in the AB0 3 perovskite films42£i 

The transition metal ions B in A2B2O7 reside on a 
three-dimensional pyrochlorc lattice [sec Fig.Qx] in which 
each B atom is centered in an octahedral cage of oxygen 
atoms. The transition metal ion B is thus subjected to 
a cubic crystalline field which splits the d-orbitals into 
a lower-lying t 2g manifold and a higher lying (typically 
on the order of a few electron volts) e g manifoldj 43 ' 81 If 
spin-orbit coupling is also present as we assume here, it 
further splits the t 2g manifold into a lower- lying j = 3/2 
manifold and higher- lying j = 1/2 manifold, with the 
size of the splitting given by the size of the spin-orbit 
coupling^ To the lowest order, the spin-orbit coupling 
does not affect the e g manifold, though at second-order 
virtual transitions to the t 2g manifold do lead to an ef- 
fective spin-orbit coupling within the e g manifold.— 
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FIG. 3. (color online) The band structure of Eq. |T} along the 
high symmetry directions [see Fig. lib] of a single kagome layer 
with t p = —2t 3 /3, t 3 = — t, t and A = 2t,4t as shown within 
the figure. Green (light gray) lines indicate filling fractions of 
of a Z2 TI. The thickness of the line indicates the size of the 
gap. In the cases shown Z2 TI occur at t2 S filling fractions: 

( a ) 9' 9' (b) 91 91 ( C ) 91 (°-) 95 gi g- 



In this work, we will focus our attention primarily on 
a transition metal ion ci-shelling filling of 6 electrons or 
less, which means we will focus on the t 2g manifold. In 
the captions of Figs.[5][()]we have used t 2g filling fractions 
for convenience of counting the bands. This must be 
converted to a total ci-shelling filling when considering 
actual materials by simply writing the t 2g filling in the 
form y/6 and then dividing y by 10. For example, a 
t 2g filling of 5/6 is a total d-shell filling of 5/10=1/2, 
and a t 2g filling of 7/9=(2*7/3)/6 is a total d-shell filling 
of (2*7/3)/10=7/15. We note that a half-filled d-shell 
corresponds to a half-filled j = 1/2 shell when spin-orbit 
coupling is present. 

Throughout this work, we will focus on 4c? and 5d tran- 
sition metal ions. For these ions, the spin orbit-coupling 
is largc-up to 2-4 times the value of the nearest neighbor 
hopping in the structures we consider . 41 i 49 The starting 
point of our study is a tight binding model of the d elec- 
trons with an atomic spin-orbit coupling. Although the 
Ad and 5d orbitals are rather localized compared to s and 
p-typc orbitals, they are still rather delocalized compared 
to 3d orbitals and their direct overlap may be apprecia- 
ble for the systems of interest here. We therefore include 
i t 46 i 48 in addition to the indirect hopping via the oxygen 
orbitals42ria 




4 r m k r 6 r m k r 



FIG. 4. (color online) The band structure of Eq. fTJ along 
the high symmetry directions [see Fig. lib] of a bilayer with 
t p = — 2t s /3, t s — —t,t and A = 2t, At as shown within the 
figure. Green (light gray) lines indicate filling fractions of of a 
Z2 TI. The thickness of the line indicates the size of the gap. 
In the cases shown Z2 TI occur at t2 g filling fractions: (a) |; 

(b) |; (d) f . 

The non-interacting Hamiltonian we consider is 

Hq= ]T tiajffC^Cjp - A^li • Sj, (1) 

where the d-orbital hopping takes the for m 46 ' 48 

tia,jP = tia,jp + ^ia,jp- (^) 

The hopping parameter Eq. ([2]) contains both an indi- 
rect and a direct hopping term between the <i-orbitals. 
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The details of the indirect hopping on the pyrochlores 
A 2 B 2 07, which is mediated by intermediate oxygen p- 
states, are given in Refs. pfJ - l42j and the details of the 
direct hopping are summarized in App. [A] The operator 
c ia creates a d-electron on site i with spin and ti g orbital 
state a, while the operator Cjp annihilates a <i-orbital 
electron on site j with spin and ti g orbital state /3. The 
matrix elements Eq. ([2]) are diagonal in the spin degree of 
freedom. In the presence of spin-orbit coupling in Eq. ([TJ , 
it is however convenient to rotate the spin quantization 
axis to the local frame, in which case the effective hop- 
ping parameters become spin-dependent. Here, A > is 
the intrinsic spin-orbit coupling in the system which acts 
within the ti g manifold so 1 = 1, and is the spin of the 
electron in a t^g <i-orbital on site i42r— In the 5d oxides, 
the strength of the spin-orbit coupling is estimated to be 
0.2-0.7 eV and the hopping strength is on the order of 
0.4-0.6 eV i 39 i 42 In order to cover the physical ranges and 
emphasize the physics associated with strong spin-orbit 
coupling, we will consider X/t in the range of 2-4, with t 
the strength of the hopping between d-orbitals mediated 
by the oxygen (indirect hopping) 42r— The hopping t will 
set our reference energy scale in this paper. 
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FIG. 5. (color online) The band structure of Eq.fTJ) along the 
high symmetry directions [see Fig |lb| of a triangular-kagome- 
triangular layer with t v — — 2t s /3, t s = —t,t and A = 2t, At 
as shown within the figure. Green (light gray) lines indicate 
filling fractions of of a Z2 TI. The thickness of the line indi- 
cates the size of the gap. In the cases shown Z2 TI occur at 
t 2g filling fractions: (a) §, if, §, if; (b) \, %, f , if; (c) ±; 

(d) — i - — 

15 ' 3 ' 5 ' 15 1 

The hopping amplitude in Eq. ^ contains a direct 
d-d hopping, tf^jp, in addition to the indirect hopping 
via the oxygen orbitalsi 46 i 48 The direct hopping is pa- 
rameterized by the strength of the c-bonds, t^, and the 
7r-bonds, t p (see App. [A"]) . Following Refs. and |48|. 
we consider a set of representative ratios to explore a 
realistic parameter space: We set t p = —2t s /3 and con- 
sider the cases of t s = —t and t s = t. For the bilayer 
system shown in Fig. [3J the phase diagrams in Fig. [7] il- 
lustrate the evolution of the ground state of the system 
as a function of t p and t s for two representative values of 



spin-orbit coupling at a total d-shcll filling of 5 electrons. 




FIG. 6. (color online) The band structure of Eq.JTJ along the 
high symmetry directions [see Fie llbj of a kagome-triangular- 
kagome layer with t p = — 2t s /3, t s — —t,t and A = 2t,At as 
shown within the figure. Green (light gray) lines indicate fill- 
ing fractions of of a Z2 TI. The thickness of the line indicates 
the size of the gap. In the cases shown Z2 TI occur at t2 g 
filling fractions: (a) g; (b) f , g; (c) if; (d) if, if. 

Our main results for the band structure of Eq. (UJ are 
shown in Figs.[3][ni Light green (gray) lines indicate filling 
fractions for which the system is a Z2 TI, and the width 
of the line is proportional to the size of the gap at that 
filling. We note that the bilayer system is not inversion 
symmetric so one must use a more general formulation^ 
to compute the Z2 invariant than the parity eigenvalues 
at time-reversal invariant momenta^ The numerical val- 
ues of the ti g filling are given in the figure captions, and 
can be converted to the total d-shell filling as described a 
few paragraphs above. This information can be used by 
experimental groups to help identify candidate materials 
for the A2B'207/A2B207/A 2 B'2C>7 sandwich structure. 
We note that in order to obtain some of the desired filling 
fractions for topological phases, it may be necessary to 
use two different A-sitc elements to create structures of 
the form k 2 W 2 Q>i I k 2{ Y_ x) k\ x ^, 2 0i I 2^ and tune 
x appropriately. This is necessary in order to obtain a 
naive "fractional" filling of the £2g-orbitals of the type 
7/9, 8/9, 3/4, etc. in which the fraction is not of the 
form n/6 for some integer n. 

A few general patterns emerge from Figs. [SUB] First, 
while for the case of t s = —t the j = 1/2 and j = 3/2 
manifolds are already clearly separated for A = 2t, the 
A — > 00 limit is well approximated by A = At because 
the separation of the upper j = 1/2 and lower j = 3/2 
manifold is large compared to the typical band widths 
already for A = At. Second, for t s = t the j = 1/2 
and j = 3/2 manifolds are not well separated even for 
A = At. Third, it appears that both the cases t s = —t 
and t s = t are equally favorable for finding topologi- 
cal phases, both in terms of the number of filling frac- 
tions that admit a Z2 TI and in terms of the sizes of 
the bulk gaps that the TI would possess. Fourth, we 
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s 

FIG. 7. (color online) Phase diagrams for the bilayer case 
in Fig. [2] at a ti g filling of 5/6 and a total d-shell filling of 
1/2 for (a) A = 2t and (b) A = At. The different colors 
represent different phases: TM is a topological metal (finite 
direct band gap, negative indirect gap, with non-trivial Zi 
index for fully occupied bands), M is a metal, I is a trivial 
insulator, and TI is a topological insulator. The black line is 
t p — —2t s /3. The * represent (-1, 2/3) and the x represent 
(1, -2/3), corresponding to the two values of t s = ±1 wc 
considered in Figs. 1 3161 



find that the simpler sandwich structure with only a sin- 
gle A-site element will only realize a Zi TI in the bi- 
layer at tig filling fraction 5/6 (a half-filled d shell, and 
a half- filled j = 1/2 manifold in the A — > 00 limit). 
This is because the bilayer is the only system with an 
even number of lattice sites in the unit cell and hence, 
the band structure can be gapped for a half-filled d 
shell. If one expands the class of systems considered to 
include the AaB^Oy/Aafi^A^BaOr/AaB'aOy-type 
structure, our results suggest that one might be able to 
find Z 2 TI with a larger bulk gap. We caution, how- 
ever, that a more accurate estimate of the bandstructure 
should be obtained within density functional theory. 

As the bilayer shown in Fig. [2] is the only system that 
exhibits a Z2 TI for for a sandwich structure of the form 
A2B'207/A2B20y/A2B'207 in the parameter ranges of 
the Hamiltonian (|T|) wc considered, we have presented 



more complete phase diagrams at tig filling 5/6 (total d 
shell and j = 1/2 manifold half-filled) in Fig. [7] Overall, 
the phase structure as a function of the parameters t s and 
t p is complicated, as might have been anticipated from 
earlier work^ However, it is clear that a larger spin-orbit 
coupling does favor expanded regions of a TI phase. 

Thus far, we have only considered the non-interacting 
Hamiltonian (Q}. We will now add interactions to the 
Hamiltonian and see how this changes the overall picture 
presented above. Our main result is that interactions will 
drive magnetic order in the system which can give rise 
to a zero magnetic field quantum Hall state known as a 
Chcrn insulator or quantum anomalous Hall state. 

III. THE INTERACTING CASES 

Many of the bulk pyrochlorc transition metal oxides ex- 
hibit magnetic phases^ and these phases are ultimately 
driven by electronic interactions. It is thus important to 
consider the effect of interactions on the non-interacting 
results we presented in Sec. HU To do so, we will add a 
Hubbard U term to Eq. JT]) and consider only the j = 1/2 
manifold by assuming the j = 3/2 manifold is well sepa- 
rated in energy (i.e., spacing is large compared to U), 

Hu = Uy^nifTin,. (3) 

where f,i refer the j z projections in the j = 1/2 mani- 
fold. In our calculations, we further assume a half-filled 
d-shell (equivalent to a half-filled j = 1/2 manifold). Wc 
study the problem within the Hartree-Fock approxima- 
tion at zero temperature. In this approximation, Hu — > 

-t / E 4 ( 2 (j' i )-j l - (j' i > 2 ) J where ji = \ E a /3=t,.L 4a CT ^ c ^- 
The ji is proportional to the spontaneous local moment 
of the d-orbitals4& We mesh the two-dimensional Bril- 
louin zone with a 150 by 150 grid of points and compute 
the band structure and the unrestricted local magnetic 
moments (j) self-consistently. Details of the calculation 
are given in App. [B] The resulting phase diagrams for a 
half-filled e?-shell are shown in Fig. [8] For gapped phases 
with zero magnetic moments (where time-reversal sym- 
metry is preserved) we compute the Z2 invariant from 
the formulation of Ref. [H| ■ For the gapped phases that 
break time-reversal symmetry, we compute the Chern 
number^ When the Chern number is nonzero (due to 
the non-trivial topological feature of the bulk bands) the 
system is guaranteed to possess gaplcss edge modes which 
support a quantum Hall effect. 

A few features of Fig. [5] are worth calling attention 
to. First, for all systems and tight-binding parame- 
ter regimes considered, there is finite range of U for 
which the system remains time- reversal invariant, while 
for U > 3.0i all systems have become magnetic, as ex- 
pected for sufficiently large U. The particular value of U 
for which this happens depends on the density of states 
at the Fermi level if the system is metallic (if larger, then 
critical U is smaller), or the magnitude of the gap if the 
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FIG. 8. (color online) The phases of the interacting model 
Q and © with (a) t s = -t, (b) t s = 0.25t, and (c) t s = t 
for a t2 g filling of 5/6 (a total d-shell and j — 1/2 filling of 
1/2) in the limit of strong spin-orbit coupling. Here M is a 
metal, I is an insulator, TI is a topological insulator, MC is a 
magnetic metal, MI is a magnetic insulator, and CI is a Chern 
insulator (zero magnetic field quantum Hall state). The most 
robust Chern insulator is in the TKT structure for t s = — t. 
Note that there is a small region of the Chern insulator in the 
bilayer system for t a — t around U/t ~ 2.5 near the end of 
the TI phase. 



system is gapped (a larger gap will have a larger critical 
U). These features can be seen by comparing the critical 
U values in Fig. [5] with the band features in Figs. [3JIB] 
As illustrated in FigJHl we observe a second-order (or 
weakly first-order) transition between the non-magnetic 
and magnetic phases in the systems that possess a Chern 
insulator phase. For the other systems (and parameter 
regimes), the magnetic transitions (not shown) most of- 
ten appear first order. 

Second, the Hubbard interaction U never has the effect 
of shifting the bands around in such a way to obtain a 
Z2 TI; if U = is not already a Z2 TI one is not later 
obtained by increasing U. Third, the bilayer is the only 
system for a half-filled d-shell that possess both a TI 




FIG. 9. (color online) Top: The magnitude of the net mag- 
netic moment, |(jt)|, and the energy gap, E g /t (E g , in d de- 
notes the indirect gap and E gt dir denotes the direct gap), as a 
function of U /t for the bilayer system with parameters corre- 
sponding to Fig. [8] Bottom: The same for the TKT system. 
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FIG. 10. (color online) The evolution of bands in a single 
kagome layer in the Hartree-Fock approximation when t 3 — t 
for (a) U = 1.90t, (b) U = 1.95t, (c) U = 2.10t and for a 
bilayer when t s = -t with (d) U = 0.65t, (e) U = 0.85t, (f) 
U = 1.05t. For the single kagome layer this corresponds to 
the M— >MC— >-MI transition in Fig. [8] and for the bilayer this 
corresponds to the I— !>MC— >MI transition. Note the Kramer's 
degeneracy at the time-reversal invariant momenta M and V 
in (d) and the difference between FK and TK' in (e) and (f). 
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FIG. 11. (color online) The evolution of the bands when t s = t 
for a bilayer with (a) U = 2.53t, (b) U = 2.554, (c) U = 
2.564, (d) U = 2.584, and for a bilayer when t a — 0.254 with 
(e) U = 1.62*, and (f) [7 = 1.724. The figures clearly show 
that for t s = t, the phase changes from (a) a TI to (b) a 
magnetic insulator (gap is too small to be visible) to (c) a 
Chern insulator and back to (d) a magnetic insulator. The 
red (gray) horizontal line in (c) reveals the gap of the Chern 
insulator. For t 3 — 0.254, the (e) TI changes to a magnetic 
insulator (f) directly. Note the Kramer's degeneracy at the 
time-reversal invariant momenta M and T in the TI in (a) and 
(e) and the difference between FK and TK'. 



and a Chern insulator. However, the TKT system for 
t s = —t is an excellent candidate for realizing the much 
sought after zero magnetic field quantum Hall state, the 
Chern insulator ( also known as quantum anomalous Hall 
state). 

In Figs. [T0lfF2l wc detail the band structure as a func- 
tion of U for different systems and tight-binding param- 
eter values to show which band features correspond to 
various transitions in the phase diagrams in Fig. [5] We 
note that the band dispersions are generally quite dif- 
ferent from their non-interacting counterparts in Figs. [3} 
[51 Figure [TO] shows the evolution of the bands around 
U ~ 2t where the M— >MC->MI transition occurs in 
the kagome layer for t s = t, and around U ~ 0.9t 
where the I— >-MC— >-MI transition occurs in the bilayer 
for t s = —t. Band features in the bilayer t s = t tran- 
sitions TI->-MI-»CI-»MI and in the bilayer t s = 0.25t 
transition TI— >MI are shown in Fig. [TT] In the case of 
the bilayer, inversion symmetry is broken resulting in a 
dispersion from T — > K being generally different from 
the dispersion from T K', so we have plotted both. 
Finally, in Fig[TOJ we show the band evolution for the 
M->MC-^CI transition in the t s = -t for a TKT layer. 

It is interesting to supplement the information about 
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=0.07t 



















M K 



(d) 2 
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FIG. 12. (color online) The evolution of the bands when t s = 
-t for a TKT layer with (a) Metal U = 0.04i, (b) Magnetic 
conductor U = 0.074, (c) Chern insulator U = 0.404, (d) 
Chern insulator U = 0.704. From (a) to (b),(c) the gap opens 
gradually, but later the gap closes again, and the tendency 
reveals near the T-point in (d). The pattern of transitions is 
M->MC->CI. 




FIG. 13. (color online) Two kinds of magnetic configurations 
obtained within the Hartree-Fock approximation for a single 
kagome layer with (a) t s = —t, U — 2.54, (b) 4 S = 4, U = 3.04 
and a bilayer with (c) 4 S = —4, U — 3.04, (d) t s = t,U = 3.04. 
For the kagome layer, two configurations exist. One (a) is 
"rotating" and keeps the Cz symmetry, while (b) is an anti- 
ferromagnetic phase without this property. For the bilayer, 
the Cz symmetry is also retained for 4 S = —4. In the other 
phase shown in (d), the magnetic moment 4 is antiparallel to 
moment 1, and moment 3 lies in the plane 413 while moment 
2 lies in plane 412. 



the changes in the band structure as a function of in- 
creasing U with figures illustrating the magnetic config- 
urations of the local moments on the c?-orbitals when a 
magnetic transition occurs. The main results are shown 
in Figs. ITOlfTIl One striking feature is that the sign of 
t s affects the magnetic order. This is evident as well in 
closely related three-dimensional studies4£ The magnetic 
moments point in different directions on different sites 
of the unit cell which reduces the overall magnetization. 
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Yet, we observe in general a net magnetic moment, see 

Fig.m 




FIG. 14. (color online) The two kinds of magnetic con- 
figurations for a TKT layer: (a) t a = -t,U = 3.0i, (b) 
t s =t,U = 3.0i and for a KTK layer: (c) t s = -t,U = 2.5t, 
(d) t s — t,U = 3.0t. In (a), the moments 4,5 are antiparallel 
to each other, and perpendicular to the 123 plane. In (b), the 
moments 4,5 are antiparallel to each other, and parallel to 
the 123 plane. One moment in 123 will lie in the same plane 
as moment 4,5, and the other two moments will point in an 
almost opposite direction. In (d) the moment 4 is parallel to 
the 123 plane. 



good growth direction. We focused on sandwich struc- 
tures where the electronically active region is confined 
to a quasi-two dimensional few-layer system as shown in 
Fig. O The analysis of the relevant tight-binding models 
for the d-shell electrons in the presence of spin-orbit cou- 
pling revealed that the bilayer system is the most promis- 
ing candidate to realize a TI phase for a half-filled d-shcll 
(as, e.g., realized in pyrochlore iridates). For the other 
systems, TI phases occur at "fractional" fillings which 
would require one to adjust the d-shell electron occupa- 
tion by considering systems of the form A 2 (i_ a; )A'2 2 ,B207 
for the electronically active region. 

The effect of the electron-electron interaction was stud- 
ied within the Hartrce-Fock mean-field theory. We found 
that the electron-electron interaction in general stabilizes 
magnetic order with non-planar local moments^ Inter- 
estingly, we have also identified parameters at intermedi- 
ate interaction strength where a Chern insulator coexists 
with magnetic order. Within our self-consistent calcula- 
tion we find that the TKT structure is the most likely 
system (among those we considered) to realize such a 
zero-magnetic field quantum Hall state. 

While our work based on semi-realistic lattice models 
revealed several scenarios to obtain topological phases 
in layered pyrochlore oxide thin films grown along the 
[111] direction, it also suggests that the expected en- 
ergy scales (i.e. the gap) of these phases are in general 
small. Moreover, both the weakly interacting phases as 
well as the type of magnetic order at larger interactions 
crucially depend on the details of the hopping param- 
eters. In the future, it is therefore desirable to incor- 
porate band structure calculations based on the density- 
functional theory to have better estimates of the hopping 
parameters and the strength of the interaction (relative 
to the band width) needed to obtain a Chern insulator 
phase. 



A final point worth emphasizing before moving to the 
conclusions is that all the magnetic orders shown in 
Figs. H3]fl~4l are non-coplanar. This condition can be fa- 
vorable for obtaining a Chern insulator, but is not strictly 
required in a mult-band system i 50 ' 51 In fact, the physics 
of the Chern insulator we obtain in our Hartree-Fock ap- 
proximation is rather similar to the large exchange limit 
studied in Ref. [S5| . In both cases, the "spin" of the elec- 
tron is forced to adiabatically follow the local moments 
residing on the lattice. 



IV. CONCLUSIONS 

In this paper, we extended previous work of perovskite 
heterostructures^r— grown along the [111] direction to 
pyrochlore heterostructures grown along [111]. Because 
(111) and equivalent planes are natural cleavage planes 
of the pyrochlore systems, it is expected that [111] is a 
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Appendix A: Tight-binding parameters for the 
direct overlap of d-orbitals in Eq.|T} 

The tight-binding hopping amplitudes used in the 
present work were derived from the nearest-neighbor hop- 
ping amplitudes on the pyrochlore lattice, including both 
oxygen-mediated (indirect) as well as direct hopping pro- 
cesses: 



liTl _i_ tdir 
ia 7 jj3 ' ia^jj: 



(Al) 



The procedure to obtain the indirect hopping t™ ^ be- 
tween the tig electrons on the pyrochlore lattice mediated 
by the intermediate oxygens was discussed in Ref. [40| . 
Due to the extended nature of the 5d orbitals, we also 
consider the direct overlap between the d orbitals which 
yields an additional contribution to the hopping matrix 
tf^jp^S. In the following, we discuss the derivation of 



dir 



Note that the oxygen octahedra enclosing the transi- 
tion metal ions arc rotated with respect to each other at 
different sites of the unit cell. The splitting of the vac- 
uum d levels into the ti g and e g manifolds results from 
the local octahedral crystal field. Hence, the ti g orbitals 
are defined with res pec t to the local frame. Following the 
convention in Ref. j4Q|, the transformation to the local 
coordinates are 



(x,y,z)Ri, 



(A2) 



where (x, y, z) are the coordinates in the global refer- 
ence frame and (xi, yi 7 Zi) are the coordinates in the local 
frame at site i. The index i = 1, 2, 3, 4 labels the sites in 
the unit cell of the pyrochlore lattice42r— Explicitly, the 
rotation matrices are given by 



(A3) 



Ri = | 


' 2/3 
-1/3 
, 2/3 


-1/3 
2/3 
2/3 


-2/3\ 
-2/3 
1/3 J 


i?2 = 


/2/3 
-2/3 
\ 1/3 


2/3 
1/3 
-2/3 


l/3\ 

2/3 

2/3/ 


R 3 = 


/ 1/3 
2/3 
V-2/3 


-2/3 
2/3 
1/3 


2/3\ 

1/3 

2/3/ 


Ra = 


' 1/3 

-2/3 
v 2/3 


-2/3 
-2/3 
-1/3 


2/3 \ 
-1/3 
-2/3/ 



(A4) 



(A5) 



(A6) 



ct2 = arccos(2/3), 
p2 = arccos(l/V5), 

72 = 2tt - arccos(-l/V5); (A8) 

a 3 = arccos(2/3), 

p3 = arccos(2/v / 5), 

73 = arccos(2/A/5); (A9) 

Q'4 = arccos(— 2/3), 

/3 4 = 2?r - arccos(2/%/5), 

74 = 2tt - arccos(-2/V5). (A10) 

The general procedure to obtain the direct overlap is to 
expand the local tig orbitals in the c?-orbitals of the global 
frame from which the hopping amplitudes are obtained 
through the use of the Slater-Koster integrals^ 

To expand the local ti g orbitals in the global frame, 
we make use of the transformation law of spherical har- 
monics under a rotation with Euler angles (a, j3, 7)^ 

Y l m (9 i ,cf> i )= Yr'ie,^,^,^^). (All) 

m'= — I 

Here, & m ,Ja,M) = e^' a S m ,Jfi)e-^i is a (21 + 
1) x (2Z + 1) dimensional rotation matrix and the elements 
d l m , m (f3) are given in Ref. [13 with 1 = 2. To proceed, 
we introduce the standard real orbitals 



d 3z 2_ r 2 

dx 2 —y 2 



I 
1 

i 

1 

71 



%~ 2 



(Y 2 



Yi) 



Yi) 



Yo 



(A12) 
(A13) 

(A14) 

(A15) 
(A16) 



Using the Euler angles in Eqs. (|A7[) - (|A10p together with 
Eq. (|A11[) and the definition of the real orbitals, we find 
that the local ti g orbitals are expanded in the following 

way 




M 



( \yz) \ 

\zx) 

\xy) 

3z 2 — r 2 
V W - y 2 ) ) 



(A17) 



For later use, we also give the corresponding Euler angles 
in the convention of Ref. [87| : 

57T 



«1 = T , 



Pi 
71 



arccos(l/3), 
3^/4; 



(A7) 



The rotation matrices Mi are 5x3 matrices given by 



Mi = 



/ -2/9 -5/9 2/9 \ 

-5/9 -2/9 2/9 

-2/9 -2/9 5/9 

2/(3V3) 2/(3V3) 4/(3^) 

V 2/3 -2/3 / 



(A18) 
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Mo = 



/ -2/9 -2/9 

2/9 5/9 

5/9 2/9 

-4/(2^3) 2/(3^3) 

V 2/3 



5/9 
-2/9 
-2/9 
-2/(3V3) 

2/3 



(A19) 



M, = 



Ma 



( 5/9 

-2/9 
2/9 

2/(3v^) 
V -2/3 



/ 5/9 
2/9 
-2/9 

2/(3V3) 
V -2/3 



2/9 
-2/9 
5/9 

-4/(3V3) 




2/9 
2/9 
-5/9 

-4/(3V3) 




-2/9 
5/9 
-2/9 
-2/(3V3) 
-2/3 / 



(A20) 



-2/9 
-5/9 
2/9 
-2/(3V3) 
-2/3 / 



(A21) 



With the above matrices, the direct hopping between the 
tig orbitals from site j to site i is written as 



d]d. 



a/3 



(A22) 



Here, the matrices Dt account for the rotation to the local 
frame in spin space as described in Refs. |40l - l4i| . Sdd(eij) 
is the 5 x 5-matrix containing the Slater-Koster integrals 
between d-orbitals along the unit direction ey between 
site j and i and is parameterized by (ddcr), (ddn), and 
(ddS). Because the <5-bonds are typically small, we only 
keep the a and 7r-bonds and introduce the parameters 
t s = t a = (dda) 
direct hopping. 



and t p = f w = (ddir) to quantify the 



Appendix B: Some details of the Hartree-Fock 
mean-field calculation 



We studied the interacting half-filled d-shell within the 
unrestricted Hartree-Fock mean field theory. A typical 



calculation starts from twenty different randomly gener- 
ated initial magnetic configurations (j) m which define the 
initial mean-field Hamiltonians via the mean-field decou- 
pled version of the interaction Eq. together with the 
non-interacting Bloch Hamiltonian Eq. (JTJ (projected to 
the j = 1/2 manifold). 

For a given filling fraction, the occupied single-particle 
states at zero temperature are easily identified on our 
meshed Brillouin zone (150 by 150). We then obtain the 
new local magnetic moments by evaluating the expecta- 
tion values (j) out on each site. This process is iterated 
until the difference between (j) m and (j) out is less than 
a fixed tolerance, in our case 



10) 



< 10" 



(Bl) 



where n — x,y,z. We compare the results from twenty 
groups of such calculations and take the one with mini- 
mum free energy. 

Generally, we find magnetic phases for U > U c and 
the magnitudes of the local moments (j) increase with U . 
However, even for U < U c , a finite magnetic moment may 
be obtained due to finite-size effects. The finite size effect 
has been identified by varying the k-mesh of the Brillouin 
zone. We found |(j)| < 0.01 to be a robust criteria for 
which the system can be considered non-magnetic. In 
these cases, we reset the magnetic moments to zero. 

Figure [5] show the magnitudes of the net magnetic mo- 
ments as function oiU/t for the bilayer and TKT sys- 
tems, respectively, 



(B2) 



where the sum runs over the sites of the unit cell. In 
addition, Fig. |H] show the direct and indirect gaps. This 
information was used to distinguish between insulating 
and metallic phases. The topological properties were 
addressed by computing the Z-i invariant for the time- 
reversal symmetric phases and the Chern number for the 
magnetic phases. 82 ' 84 
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